Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 20 April 2009 (MN M£X style file v2.2) 



On the derivation of dynamical masses of the stellar 
clusters in the circumnuclear region of NGC 2903 

Guillermo F. Hagele 1 ' 2 *, Angeles I. Diaz 1 , Monica V. Cardaci 1 ' 3 !, Elena Terlevich 4 J 
and Roberto Terlevich 4 ^ 

1 Departamento de Fisica Teorica, C-XI, Universidad Autdnoma de Madrid, 28049 Madrid, Spain 

2 Facultad de Cs Astronoraicas y Geopsicas, Universidad Nacional de La Plata, Paseo del Bosque s/n, 1900 La Plata, Argentina 

3 XMM Science Operations Centre, European Space Astronomy Centre of ESA, P.O. Box 50727, 28080 Madrid, Spain 

4 INAOE, Tonantzintla, Apdo. Postal 51, 72000 Puebla, Mexico 



20 April 2009 



ABSTRACT 

Gas and star velocity dispersions have been derived for four circumnuclear star-forming 
regions (CNSFRs) and the nucleus of the spiral galaxy NGC 2903 using high resolution 
spectroscopy in the blue and far red. Stellar velocity dispersions have been obtained 
from the Call triplet (CaT) lines at AA 8494, 8542, 8662 A, using cross-correlation 
techniques while gas velocity dispersions have been measured by Gaussian fits to the 
H/? A 4861 A line. 

The CNSFRs, with sizes of about 100 to 150 pc in diameter, show a complex 
structure at the Hubble Space Telescope (HST) resolution, with a good number of 
subclusters with linear diameters between 3 and 8 pc. Their stellar velocity dispersions 
range from 39 to 67 kms -1 . These values, together with the sizes measured on archival 
HST images yield upper limits to the dynamical masses for the individual star clusters 
between 1.8 and 8.7 x 1O 6 M and upper limits to the masses for the whole CNSFR 
between 4.9 x 10 6 and 4.3 x 1O 7 M . 

The masses of the ionizing stellar population responsible for the Hn region gaseous 
emission have been derived from their published Ha luminosities and are found to be 
between 1.9 and 2.5 x 10 6 Mq for the star-forming regions, and 2.1 x 10 5 M Q for the 
galaxy nucleus therefore constituting between 1 and 4 per cent of the total dynamical 
mass. 

In the CNSFR star and gas velocity dispersions are found to differ by about 
20kms~ 1 with the H/3 lines being narrower than both the stellar lines and the 
[Oiii]A5007A lines. The ionized gas kinematics is complex; two different kinemat- 
ical components seem to be present as evidenced by different widths and Dopplcr 
shifts. 

The line profiles in the spectra of the galaxy nucleus, however, are consistent with 
the presence of a single component with radial velocity and velocity dispersion close 
to those measured for the stellar absorption lines. 

The presence and reach of two distinct components in the emission lines in ionized 
regions and the influence that this fact could have on the observed line ratios is of 
major interest for several reasons, among others, the classification of the activity in the 
central regions of galaxies, the inferences about the nature of the source of ionization 
for the two components and the derivation of the gaseous chemical abundances. 

Key words: HII regions - galaxies: individual: NGC 2903 - galaxies: kinematics and 
dynamics - galaxies: starburst - galaxies: star clusters. 
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1 INTRODUCTION 

The inner (~1 Kpc) parts of some spiral galaxies show high 
star formation rates and this star formation is frequently 
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arranged in a ring or pseudo-ring pattern around their nu- 
clei. This fact seems to correlate with the presence of bars 



( Roberts et al.||1979| [Friedli fc Benz|[l9 95 Saka moto et al. 
1999||Sheth et al.|2005| |Jogee et al.|2005|) and, in fact, com- 



puter models which simulate the behaviour of gas in galactic 
potentials have shown that nuclear rings may appear as a 
consequence of matter infall owing to resonances present at 
the bar edges ( Combes fc Gerin|19 85 Athanassoula 19921. 

These circumnuclear star-forming regions (CNSFRs) 
are expected to be amongst the highest metallicity regions as 
corresponds to their position near the galactic center. This 
fact makes the analysis of these regions complicated since, 
in general, their low excitation makes any temperature sen- 
sitive line too weak to be measured ( Diaz et aL||2007 \ , par- 
ticularly against a strong underlying stellar continuum. As 
was pointed out in Hagele et al. (2007; hereafter Paper I) 
the [Oiii] emission lines in CNSFRs are generally very weak, 
and in some cases unobservable. The low value of these col- 
lisionally excited lines can be explained by their over-solar 



metal abundances (e.g. Boer & Schulz 1993 1. The equivalent 
width of the emission lines are lower than those shown by 



the disc Hn regions (see for example Kennicutt et al.||1989 



Bresolin & Kennicutt 1997 Bresolin et al.|1999 1 . Combining 
GEMINI data and a grid of photo-ionization models |Dors| 
et al. ( 2008 1 conclude that the contamination of the con- 



tinua of CNSFRs by underlying contributions from both old 
bulge stars and stars formed in the ring in previous episodes 
of star formation (10-20 Myr) yield the observed low equiv- 
alent widths. 

CNSFRs, with sizes going from a few tens to a few hun- 



dreds of parsecs (e.g. Diaz et al. 2000 1 seem to be made 



of several Hn regions ionized by luminous compact stellar 
clusters whose sizes, as measured from high spatial resolu- 
tion Hubble Space Telescope (HST) images, are seen to be of 
only a few parsecs. Although these Hn regions are very lu- 
minous (M„ between -12 and -17) not much is known about 
their kinematics or dynamics for both the ionized gas and 
the stars. It could be said that the worst known parameter 
of these ionizing clusters is their mass. There are different 
methods to estimate the mass of a stellar cluster. As derived 
with the use of population synthesis models their masses 
suggest that these clusters are gravitationally bound and 
that they might evolve into globular cluster configurations 
( Maoz et al.|1996| . Classically one assumes that the system 
is virialized and determines the total mass inside a radius 
by applying the virial theorem to the observed velocity dis- 
persion of the stars (<r»). The stellar velocity dispersion is 
however hard to measure in young stellar clusters (a few 
million-years old) due to the shortage of prominent stellar 
absorption lines. The optical continuum between 3500 and 
7000 A shows very few lines since the light at these wave- 
lengths is dominated by OB stars which have weak absorp- 
tion lines at the same wavelengths of the nebular emission 
lines (Balmer H and Hel lines). As pointed out by several au- 
thors (e.g. |Ho fc Filippenko|| 1996a I, at longer wavelengths 



( 8500 A) the contamination due to nebular lines is much 
smaller and since red supergiant stars, if present, dominate 
the near-infrared (IR) light where the Can AA8498, 8542, 
8662 A triplet (CaT) lines are found, these should be easily 



observable allowing the determination of a* ( Terlevich et al. 
1990||Prada et al.|1994 i. 



as a SB(rs)bc by de Vaucouleurs et al. ( 1991 1 . Its coordinates 
are_02ooo=09' 1 32 m Vf.l, <5 2 ooo=+21° SWOT \de Vaucouleurs 
et al.|1991 l. These authors derived a Bt (mag) equal to 9.7 
and an E(B-V) ga i (mag) of 0.031. We adopt the distance 



derived by Bottinelli et al. ( 1984 1 which is equal to 8.6 Mpc, 
giving a linear scale of ~ 42 pc arcsec -1 . 

The Paschen a image obtained with the Hubble Space 
Telescope (HST) reveals in NGC 2903 the presence of a nu- 
clear ring-like morphology with an apparent diameter of ap- 
proximately 15" = 625pc ( Alonso-Herrero et al.|2 001). This 
structure is also seen, though less prominent, in the Ha ob- 



servations from Planesas et al. (1997). A large number of 



stellar clusters are identified on high resolution infrared im- 
ages in the K' and H bands, which do not coincide spatially 
with the bright Hn regions. A possible interpretation of this 
is that the stellar clusters are the result of the evolution 



of giant Hn regions (e.g. Alonso-Herrero et al. 20011. The 
global star formation rates in the nuclear ring, as derived 
from its Ha lum inosity is found to be 0.1 Mq yr -1 by 



sas et al 



(19971 and 0.7M Q yr -i by 



Plane- 



Alonso-Herrero et al. 



(20011. From CO emission observations, Planesas et al. de- 
rive a mass of molecular gas (H2) of 1.8 x 10 8 Mq inside a 
circle 1 Kpc in diameter. 

This is the second paper of a series to study the pe- 
culiar conditions of star formation in circumnuclear regions 
of early type spiral galaxies, in particular the kinematics 
of the connected stars and gas. In this paper we present 
high-resolution far-red spectra and stellar velocity dispersion 
measurements (<t») along the line of sight for four CNSFRs 
and the nucleus of the spiral galaxy NGC 2903. We have also 
measured the ionized gas velocity dispersions (<r 9 ) from high- 
resolution blue spectra using Balmer H/3 and [Oiii] emission 
lines. The comparison between <r* and a g might throw some 
light on the yet unsolved issue about the validity of the grav- 
itational hypothesis for the origin of the supersonic motions 



NGC 2903 (UGC 5079) is a well studied galaxy classified 



observed in the ionized gas in Giant Hn regions (Melnick, 
|Tenorio-Tagle fc Terlevich|p99| ). 

In Section 2 we describe the observations and data re- 
duction. We present the results in Section 3, the dynamical 
mass derivation in Section 4 and the ionizing star cluster 
properties in Section 5. We discuss all our results in Section 
6. Finally, summary and conclusions are given in Section 7. 



2 OBSERVATIONS AND DATA REDUCTION 

The data were acquired in 2000 February using the two arms 
of the ISIS spectrograph attached to the 4.2-m William Her- 
schel Telescope (WHT) of the Isaac Newton Group (ING) 
at the Roque de los Muchachos Observatory on the Spanish 
island of La Palma. The CCD detectors EEV12 and TEK4 
were used for the blue and red arms with a factor of 2 binning 
in both the "x" and "y" directions with resultant spatial res- 
olutions of 0.38 and 0.36 arcsec pixel -1 for the blue and red 
configurations respectively. The H2400B and R1200R grat- 
ings were used to cover the wavelength ranges from 4779 to 
5199 A (A c = 4989 A) in the blue and from 8363 to 8763 A 
(A c = 8563 A) in the red with resultant spectral dispersions 
of 0.21 and 0.39 A per pixel respectively, providing a com- 
parable velocity resolution of about 13 km s -1 . A slit width 
of 1 arcsec was used which, combined with the spectral dis- 
persions, yielded spectral resolutions of about 0.4 and 0.7 A 
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Table 1. Journal of Observations 



Date 


Spectral range 


Disp. 


tj a 

FWHM 


Spatial res. 


PA 


Exposure Time 


seeing FW HM 




(A) 


(Apx- 1 ) 




("px- 1 ) 


(°) 


(sec) 


(") 


2000 February 4 


4779-5199 


0.21 


~ 12500 


0.38 


50 


3 X 1200 


1.2 




8363-8763 


0.39 


~ 12200 


0.36 


50 


3 x 1200 




2000 February 5 


4779-5199 


0.21 


~ 12500 


0.38 


345 


3 x 1200 


1.6 




8363-8763 


0.39 


~ 12200 


0.36 


345 


3 x 1200 





°%WHM — V^FWHM 



FWHM in the blue and the red, respectively, measured on 
the sky lines. Table [I] summarizes the instrumental configu- 
ration and observation details. Several bias and sky flat field 
frames were taken at the beginning and end of each night in 
both arms. In addition, two lamp flat field and one calibra- 
tion lamp exposures per telescope position were performed. 
The calibration lamp used was CuNe+CuAr. 

The data were processed and analysed using IRAE^Jrou- 
tines including the usual procedures of removal of cosmic 
rays, bias subtraction, division by a normalized flat field 
and wavelength calibration. Further details concerning each 
of these steps can be found in Paper I. 

Neither atmospheric extinction nor flux calibration were 
performed, since our purpose was to measure radial veloci- 
ties and velocity dispersions. With this purpose, spectra of 
11 template velocity stars were acquired to provide good 
stellar reference frames in the same system as the galaxy 
spectra for the kinematic analysis in the far-red. They cor- 
respond to late-type giant and supergiant stars which have 
strong CaT features (see Diaz et al. 19891. The spectral 



types, luminosity classes and dates of observation of the stel- 
lar reference frames used as templates are listed in Table 2 
of Paper I. 



3 RESULTS 

Two different slit positions (SI and S2) were chosen in or- 
der to observe the nucleus of the galaxy and 4 CNSFRs. 
Fig. [l] shows them, superimposed on photometrically cali- 
brated optical and infrared images of the circumnuclear re- 
gion of NGC 2903 acquired with the Wide Field and Plane- 
tary Camera 2 (WFPC2; PCI) and the Near-Infrared Cam- 
era and Multi-Object Spectrometer (NICMOS) Camera 3 
(NIC3) on board the HST. These images have been down- 
loaded from the Multimission Archive at STScI (MASTQ 
The optical image was obtained through the F606W (wide 
V) filter, and the near-IR one, through the F160W (H). The 
CNSFRs have been labelled following the same nomencla- 
ture as in |Planesas et allj ( |1997[ ) from their identifications on 
Ha maps. The plus symbol in both panels of Fig. [I] repre- 
sents the position of the nucleus as given by the Two Micron 
All S ky Survey Team, 2M ASS Extended Ojects - Final Re- 
lease ( Jarrett et al.|2003 |. 

Fig. [2j shows the spatial profiles in the H/3 and 



1 IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Sci- 
ence Foundation (NSF). 

2 http:/ /archive. stsci.edu/hst/wfpc2 



[Om] 5007 A emission lines (upper and middle panels) and 
the far-red continuum (lower panels) along each slit position. 
The emission line profiles have been generated by collapsing 
11 pixels of the spectra in the direction of the resolution at 
the central position of the lines in the rest frame, A 4861 and 
A 5007A respectively, and are plotted as dashed lines. Con- 
tinuum profiles were generated by collapsing 11 resolution 
pixels centred at 11 A to the blue of each emission line and 
are plotted as dash-dotted lines. The difference between the 
two, shown by solid lines, corresponds to the pure emission. 
The far-red continuum has been generated by collapsing 11 
pixels centred at A8620A. 

The regions of the frames to be extracted into one- 
dimensional spectra corresponding to each of the identified 
CNSFRs, were selected on these profiles with reference to 
the continuum emission both in the blue and in the red. 
These regions are marked by horizontal lines and labelled in 
the figure. In the H/3 profile of slit position SI, an almost 
pure emission knot between regions R1+R2 and R7 can be 
seen with a very weak continuum. We labelled this region as 
X. In the case of the [Om] 5007 A emission line profile, the 
main contributor comes from the continuum since the actual 
emission line is very weak as expected for high metallicity 
Hn regions ( Diaz et al.|2007 l. Spectra in slit position SI are 
placed along the circumnuclear regions located to the North 
and North-East of the nucleus and therefore any contribu- 
tion from the underlying galaxy bulge is difficult to assess. In 
the case of slit position S2 which crosses the galactic nucleus, 
this fact can be used to estimate that contribution. For the 
blue spectra the light from the underlying bulge is almost 
negligible amounting to, at most, 5 per cent at the H/3 line. 
For the red spectra, the bulge contribution is more impor- 
tant. From Gaussian fits to the A 8620 A continuum profile 
we find that it amounts to about 15 per cent for the lowest 
surface brightness region, R1+R2. On the other hand, the 
analysis of the broad near-IR HST-NICMOS image shown 
in Fig. [l] indicates less contrast between the emission from 
the regions and the underlying bulge which is very close to 
the image background emission. Its contribution is about 25 
per cent for the weak regions, R2 and R7 along position slit 
SI. 

Figs. [3] and [4] show the spectra of the observed regions 
and the nucleus, respectively, split into two panels corre- 
sponding to the blue and the red spectral range. The blue 
spectra show the Balmer H/3 recombination line and the 
weak collisionally excited [Om] lines at AA4959, 5007 A. The 
red spectra show the stellar Call triplet lines in absorp- 
tion (CaT) at AA 8498, 8542, 8662 A. In some cases, as for 
example R1+R2, these lines are contaminated by Paschen 
emission which occurs at wavelengths very close to those of 
the CaT lines. Other emission features, such as OlA8446, 
[Cln] A 8579, Pa 14 and [Fen] A 8617 are also present. In this 
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Figure 1. Left: F606W (wide V) image centred on NGC 2903 obtained with the WFPC2 camera (PCI) of the HST. Right: HST-NICMOS 
(NIC3) image obtained through the F160W filter. For both images the orientation is north up, east to the left. The location and P.A. of 
the WHT-ISIS slit positions, together with identifications of the CNSFRs extracted, are marked. 



case, a single Gaussian fit to the emission lines was per- 
formed and the lines were subsequently subtracted after 
checking that the theoretically expected ratio between the 
Paschen lines was satisfied. The observed red spectra of 
R1+R2 are plotted in Fig. [3] with a dashed line. The solid 
lines show the subtracted spectra. 

Fig. [5] shows the spectrum of the almost pure emission 
knot marked in the profiles of NGC 2903, region X. In this 
case the blue range of the spectrum presents very intense 
emission lines, while in the red range we can appreciate a 
spectrum similar to those shown by the CNSFRs studied 
by Planesas et al. (1997) suggesting a region with similar 
characteristics but with low continuum surface brightness. 



3.1 Kinematics of stars and ionized gas 

A detailed description of the methods and techniques used to 
derive the values of radial velocities and velocity dispersions 
as well as sizes, masses and emission line fluxes has been 
given in Paper I. Therefore only a brief summary is given 
below. 



Stellar analysis 

Stellar radial velocities and velocity dispersions were ob- 
tained from the conspicuous CaT absorption lines using the 
cross-correlation technique, described in detail by |Tonry fc| 
|Davis| ( |1979[ ) . This method requires the comparison with a 
stellar template that represents the stellar population that 
best reproduces the absorption features. This has been built 
from a set of 11 late- type giant and supergiant stars with 
strong CaT absorption lines. We have followed the work by 



Nelson & Whittle ( 1995 1 with the variation introduced by 
Palacios et al.| ( |1997| ) of using the individual stellar tem- 
plates instead of an average. This procedure will allow us to 



correct for the known possible mismatches between template 
stars and the region composite spectrum. The implementa- 
tion of the method in the external package of IRAF XCSAO 
(Kurtz & Mink 19981 has been used, as explained in Paper 
I. 

To determine the line-of-sight stellar velocity and veloc- 
ity dispersion along each slit, extractions were made every 
two pixels for slit position SI and every three pixels for slit 
position S2, with one pixel overlap between consecutive ex- 
tractions in this latter case. In this way the S/N ratio and 
the spatial resolution were optimised. The stellar velocity 
dispersion was estimated at the position of each CNSFR 
and the nucleus using an aperture of five pixels in all cases, 
which corresponds to 1.0 x 1.8arcsec 2 . The velocity disper- 
sion (a) of the stars (a*) is taken as the average of the a 
values found for each stellar template, and its error is taken 
as the dispersion of the individual values of a and the rms of 
the residuals of the wavelength fit. These values are listed in 
column 3 of Table [2] along with their corresponding errors. 

The radial velocities have been determined directly from 
the position of the main peak of the cross-correlation of each 
galaxy spectrum with each template in the rest frame. The 
average of these values is the final adopted radial velocity. 



Ionized gas analysis 

The velocity dispersion of the ionized gas was estimated for 
each observed CNSFR and for the galaxy nucleus from gaus- 
sian fits to the H/3 and [Oiii] A 5007 A emission lines using 
five pixel apertures, corresponding to 1.0 x 1.9arcsec 2 . For 
a single Gaussian fit, the position and width of a given emis- 
sion line is taken as the average of the fitted Gaussians to the 
whole line using three different suitable continua ( |Jimenez^| 
Benito et al.|2000 \ , and their errors are given by the disper- 



Dynamical masses of the CNSF clusters of NGC 2903 5 




si 




S2 





Figure 2. Spatial profiles of H/3, [Om] A5007A and the far red (upper, middle and lower panels respectively) for each slit. For the 
emission lines, the profiles correspond to line+continuum (dashed line), continuum (dashed-dotted line) and the difference between them 
(solid line), representing the pure emission from H/3 and [Om] respectively. For the far red profiles, the solid lines represent the continuum. 
Pixel number increases to the North. Horizontal small lines show the location of the CNSFRs and nuclear apertures. 
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Figure 3. Blue (left) and red (right) rest frame normalized spectra of the observed CNSFRs. For R1+R2, the dashed 
obtained spectrum; the solid line represents the spectrum after subtracting the emission lines (see text). Notice the absence 
emission lines in the red spectral range for R4 and R7. 
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Figure 4. Blue (left) and red (right) rest frame normalized spectra of the nucleus. 
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Figure 5. Idem as Fig. [3] for region X. The dashed line shows the obtained spectrum; the solid line represents the spectrum after 
subtracting the emission lines (see text). 



sion of these measurements taking into account the rms of 
the wavelength calibration. 

In all the studied regions, however, the best fit for the 
H/3 line is obtained with two different components having 
different radial velocities of up to 35 km s _1 . We used the 
widths of those components as an initial approximation to 
fit the [Om] lines which, due to their intrinsic weakness, 
show lower S/N ratio, and found them to provide also an 
optimal two component fit. The radial velocities found for 
the narrow (broad) components of both H/3 and [Om] are 
the same within the errors. An example of the two-Gaussian 
fit is shown in Fig. [6] 

For each CNSFR, the gas velocity dispersions for the 
H/3 and [Om] A 5007 A lines derived using single and double 
line Gaussian fit, and their corresponding errors are listed in 
Table [2] Columns 4 and 5, labelled 'One component', give 
the results for the single Gaussian fit. Columns 6 and 7, 
and 8 and 9, labelled 'Two components - Narrow' and 'Two 
components - Broad', respectively, list the results for the two 
component fits. The last column of the table, labelled Av„t, 
gives the velocity difference between the narrow and broad 
components. This is calculated as the average of the H/3 and 
[Om] fit differences. Taking into account the errors in the 
two component fits, the errors in these velocity differences 
vary from 10 to 15kms _1 . 

We have also determined the distribution along each 



slit position of the radial velocities and the velocity disper- 
sions of the ionized gas using the same procedure as in the 
case of the stars, that is using spectra extracted every two 
pixels for SI and every three pixels, superposing one pixel 
for consecutive extractions, for S2. These spectra, however, 
do not have the required S/N ratio to allow an acceptable 
two-component fit, therefore a single-Gaussian component 
has been used. The goodness of this procedure is discussed 
below in Section 4. 



3.2 Emission line ratios 

We have used two different ways to integrate the inten- 
sity of a given line: (1) if an adequate fit was attained by 
a single Gaussian, the emission line intensities were mea- 
sured using the SPLOT task in IRAF. For the H/3 emis- 
sion lines a conspicuous underlying stellar population is 
inferred from the presence of absorption features that de- 
press the lines (e.g. see discussion in Diaz 1988). Exam- 
ples of this effect can be appreciated in Fig. |7| We have 
defined a pseudo-continuum at the base of the line to 
measure the line intensities and minimize the errors intro- 



duced by the underlying population (for details see Hiigele 
et al. 2006 1. (2) When the optimal fit was obtained by 
two Gaussians the individual intensities of the narrow and 
broad components are estimated from the fitting parameters 
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Table 2. Velocity dispersions. 
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Figure 6. Sections of the normalized spectrum of R1+R2 (solid line). The left panel shows from 4861 to 4876 A containing H/3 and 
the right panel shows from 5007 to 5022 A containing the [Oiii]A5007A emission line. For both we have superposed the fits from the 
ngaussfit task in IRAF; the dashed-dotted line is the broad component, the dotted line is the narrow component and the dashed line is 
the sum of both. 



(I = 1.0645 A x FWHM = V2tt A a; where I is the Gaussian 
intensity, A is the amplitude of the Gaussian, FWHM is the 
full width at half-maximum and a is the dispersion of the 
Gaussian). A pseudo-continuum for the H/3 emission line was 
also defined in these cases. The statistical errors associated 
with the observed emission fluxes have been calculated with 
the expression <t ; =ct c N 1/2 [1 + EW/(NA)] 1/2 , where o~i is 
the error in the observed line flux, a c represents the stan- 
dard deviation in a box near the measured emission line and 
stands for the error in the continuum placement, N is the 
number of pixels used in the measurement of the line inten- 
sity, EW is the line equi valent width and A is the wavele ngth 
dispersion in A pixel -1 ( Gonzalez- Delgado et al.|l994 k For 
the H/3 emission line we have doubled the derived error, ai, 
in order to take into account the uncertainties introduced 



by the presence of the underlying stellar population ( Hagele 
et al.|2006| >. 



The logarithmic ratio between the emission line intensi- 
ties of [Oiii] A 5007 A and H/3 and their corresponding errors 
are presented in Table [3] We have also listed the logarith- 
mic ratio between the emission line fluxes of [No] A 6584 A 
and Hcv together with their corresponding e rrors from|Diaz 



et al. 



(20071 for R1+R2 and R4, and from Planesas et al 



( |1997 1 for R7 and the nucleus. For these last two objects, the 
[Nil] and Hcv emission are derived from the Ha+[Nll] nar- 
row band images. Planesas and collaborators estimated the 
relative contribution of the [Nil] AA 6548,6584 A lines from 
the Ha and [Nil] A 6584 A equivalent width measurements of 



Stauffer| fi"982j ). The logarithmic ratios for R1+R2 and R4 



derived using this procedure are similar (-0.64 ±0.04 and - 
0.65 ±0.04). For the X region we assumed a value of -0.38 
(without error) for this ratio as given by Diaz et al. (2007) 
for the other CNSFRs of this galaxy. 



4 DYNAMICAL MASS DERIVATION 

The mass of a virialized stellar system is given by two pa- 
rameters: its velocity dispersion and its size. 

In order to determine the sizes of the stellar clusters 
within our observed CNSFRs, we have used the retrieved 
wide V HST image which provides a spatial resolution of 
0.045 "per pixel. At the adopted distance for NGC 2903, this 
resolution corresponds to 1.9pc/pixel. Fig.|8]shows enlarge- 
ments around the different regions studied with intensity 
contours overlapped. 

We find, as expected, that the four CNSFRs studied 
here are formed by a large number of individual star-forming 
clusters. Regions Rl and R2 form an interrelated complex 
with two main structures: Rl and R2 surrounded by, at least, 
31 secondary structures. Rl is resolved into two separate 
clusters, labelled Rla and Rib and shows other two associ- 
ated knots, labelled Rlc and Rid. R2 is made up of a main 
cluster, labelled R2, and four secondary ones labelled "a" to 
"d" . The rest of the knots can not be directly related to any 
of the two main structures, and have been labelled as R12 
followed by a letter from "a" to "v" . The same criterion has 
been applied for regions R4 and R7 for which we find 28 and 
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Figure 7. Enlargement of the blue rest frame normalized spectra of R1+R2 (left) and R7 (right). 



Table 3. Line ratios. 







One component 


Two components 










Narrow 


Broad 




Region 


Slit 


log([Om]5007/Hy8) 


log([Om]5007/H/3) 


log([Oiii]5007/H/3) 


log([Nn]6584/Ha) 


R1+R2 


SI 


-1.03±0.05 


-1.64±0.Q9 


-0.76±0.10 


-0.37±0.01 a 


R1+R2 


S2 


-0.92±0.07 


-1.51±0.11 


-0.57±0.14 




R4 


S2 


-0.78±0.11 


-I.01i0.10 


-0.71±0.15 


-0.38±0.01 a 


K7 


SI 


-0.78±0.10 


-1.79±0.12 


-0.40±0.18 


-0.66±0.04 b 


N 


S2 


-0.50±0.07 


-0.92±0.10 


-0.23±0.18 


-0.68±0.04 6 


X 


SI 


-0.93±0.06 


-1.45±0.Q9 


-0.77±0.11 


-0.38: c 



a from Diaz et al. 1 2007 1 
6 from Vlanesas et al. (1997) 

c Assumed from the values given by Diaz et al. (2007) for the other observed CNSFR in NGC 2903 



27 individual clusters, respectively. On the other hand, the 
galaxy nucleus seems to be formed by one main structure 
with two small secondary ones to the South West, and some 
very weak knots which are difficult to appreciate in the pic- 
ture. The secondary structures have been designated by Na 
and Nb. In all cases the knots have been found with a de- 
tection level of 10 a above the background. All these knots 
are within the radius of the regions defined by Planesas et 
al. (1997). We have to remark that our search for knots has 
not been exhaustive since that is not the aim of this work. 

We have fitted circular regions to the intensity contours 
corresponding to the half light brightness distribution of 
each single structure (see Fig. [8|, following the procedure 
given in |Meurer et al.| ( |1995[ ), assuming that the regions 
have a circularly symmetric Gaussian profile. The radii of 
the single knots vary between 1.5 and 4.0 pc. Table [4] gives, 
for each identified knot, the position, as given by the as- 
trometric calibration of the HST image; the radius of the 
circular region defined as described above together with its 
error; and the peak intensity in counts, as measured from 
the WFPC2 image. The nucleus of NGC 2903 is rather com- 
pact, but resolved, with a radius of 3.8 pc. Obviously, the 
spatial resolution and, at least, the estimates of the smaller 
radii depend on the distance to each galaxy. 

Upper limits to the dynamical masses (M*) inside the 
half-light radius (R) for each observed knot have been 
estimated under the following assumptions: (i) the sys- 
tems are spherically symmetric; (ii) they are gravitationally 
bound; and (iii) they have isotropic velocity distributions 
[a 2 (total) = 3 a 2 ]. The general expression for the virial mass 



of a cluster is r\ a 2 R/G, where R is the effective gravitational 
radius and rj is a dimensionless number that takes into ac- 
count departures from isotropy in the velocity distribution 
and the spatial mass distribution, binary fraction, mean sur- 
face d ensity, etc. (|Boily et al.|[2005] [Fleck et al.||2006| ). Fol- 



lowing Ho Sz Filippenko| (|1996a|b[), and for consistence with 



Paper I and Hagele| ( |2008[ ), we obtain the dynamical masses 
inside the half-light radius using r\ = 3 and adopting the half- 
light radius as a reasonable approximation of the effective 
radius. Other authors (e.g. Spitzer|1987 Smith & Gallagher 
|2001| |Moll et aL"1|2008[ ) assumed that the r\ value is about 
9.75 thus obtaining the total mass. On the absence of any 
knowledge about the tidal radius of the clusters, we adopted 
this conservative approach. On the derived masses, the dif- 
ferent adopted rj values act as multiplicative factors. 

It must be noted that we have measurements for the 
size of each knot, but we do not have direct access to 
the stellar velocity dispersion of each individual cluster, 
since our spectroscopic measurements encompass a wider 
area (1.0 x 1.9arcsec 2 , which corresponds approximately to 
42 x 76 pc 2 at the adopted distance of NGC 2903) that in- 
cludes the whole CNSFRs to which each group of knots be- 
long. 

The estimated dynamical masses for each knot and their 
corresponding errors are listed in Table [5] For the regions 
that have been observed in more than one slit position, we 
list the derived values using the two separate stellar veloc- 
ity dispersions. The dynamical masses in the rows labelled 
"sum" have been found by adding the individual masses in 
a given CNSFR: R12, R4 and R7, as well as the galaxy 
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nucleus, N. We have taken the average of the R12sum esti- 
mated from slit positions SI and S2 as the "adopted" dy- 
namical mass for R12. The fractional errors of the dynamical 
masses of the individual knots and of the CNSFRs are listed 
in column 4. 



5 IONIZING STAR CLUSTER PROPERTIES 

For each of the CNSFR: R12, R4 and R7, the total number 
of ionizing photons was derived from the total observed Ha 
luminosities given by Planesas et al. (1997), correcting for 



Dynamical masses of the CNSF clusters of NGC 2903 11 

Table 4. Positions, radii and peak intensities derived from the HST-F606W image. 



Region 


position 


R 


I 


Region 


position 


R 


I 




a J2000.0 


<5 72000.0 


(pc) 


(counts) 




Q ./2000.0 


<5j2000.0 


(PC) 


(counts) 


Rla 


9 11 32 m lffi06 


+21°30W86 


1.6±0.1 


4083 


R4n 


9 h 32 m 10?09 


+21°29'59'.'86 


1.5±0.3 


240 


Rib 


9 h 32 m 10?06 


+21°30W95 


1.6±0.2 


4015 


R4o 


9 h 32 m 10!18 


+21°29'59'.'70 


2.3±0.3 


237 


Rlc 


9 h 32 m 10?07 


+21°30'07('10 


1.5±0.1 


363 


R4p 


9 h 32 m 10?16 


+21°29'59'.'10 


1.6±0.2 


236 


Rid 


9 h 32 m 10f05 


+21°30W08 


2.1±0.1 


254 


R4q 


9 h 32 m 10!25 


+21°29'58'.'92 


2.7±0.3 


228 


R2 


9 h 32 m 09?98 


+21°30W65 


1.9±0.1 


2892 


R4r 


9 h 32 m 10?10 


+21°29'59'.'71 


1.7±0.2 


227 


R2a 


9 h 32 m 09?96 


+21°30'06?28 


1.9±0.1 


390 


R4s 


9 h 32 m 10!21 


+21°30'00'.'18 


2.7±0.3 


216 


R2b 


9 h 32 m 09?96 


+21°30W61 


1.9±0.2 


307 


R4t 


9 h 32 m 10?08 


+21°29'59'.'83 


1.8±0.2 


208 


R2c 


9 h 32 m 09?96 


+21°30W00 


1.9±0.1 


209 


R4u 


9 h 32 m 10!ll 


+21°29'58'.'94 


2.4±0.3 


198 


R2d 


9 h 32 m 09?98 


+21°30W30 


1.8±0.2 


206 


R4v 


9 h 32 m 10?20 


+21°30'00'.'37 


2.1±0.3 


196 


R12a 


9 h 32 m 09?96 


+21°30'07y68 


1.9±0.1 


995 


R4w 


9 h 32 m 10?23 


+21°30'00'.'67 


2.0±0.2 


195 


R12b 


9 h 32 m 10?09 


+21°30'05?94 


2.5±0.1 


580 


R4x 


9 h 32 m 10!22 


+21°30'01'.'17 


3.2±0.4 


194 


R12c 


9 h 32 m 10 s 08 


+21°30W03 


2.5±0.1 


574 


R4y 


9 h 32 m 10!12 


+21°29'58'.'82 


1.9±0.3 


187 


R12d 


9 h 32 m 10f08 


+21°30W42 


2.1±0.1 


427 


R4z 


9 h 32 m 10!19 


+21°30'00'.'25 


2.0±0.3 


180 


R12e 


9 h 32 m 10 s 15 


+21°30W80 


1.6±0.1 


381 


R4o 


9 h 32 m 10!24 


+21°30'00'.'43 


2.0±0.3 


173 


R12f 


9 h 32 m 10?03 


+21°30W38 


1.8±0.1 


341 


R7 


9 h 32 m 10?38 


+21°30'10'.'93 


2.7±0.1 


984 


R12g 


9 h 32 m 10?04 


+21°30W51 


3.2±0.2 


293 


R7a 


9 h 32 m 10?32 


+21°30'09'.'74 


2.7±0.2 


715 


R12h 


9 h 32 m 10?04 


+21°30W91 


3.6±0.4 


266 


R7b 


9 h 32 m 10!27 


+21°30'09'.'99 


3.4±0.1 


579 


R12i 


9 h 32 m 10 s 15 


+21°30W45 


3.2±0.3 


255 


R7c 


9 h 32 m 10?34 


+21°30'09'.'65 


2.7±0.3 


520 


R12j 


9 h 32 m 10 s 13 


+21°30W42 


3.2±0.4 


238 


R7d 


9 h 32 m 10?39 


+21°30'09'.'17 


2.6±0.2 


470 


R12k 


9 h 32 m 10?ll 


+21°30W54 


4.0±0.5 


233 


R7e 


9 h 32 m 10?26 


+21°30'09'.'56 


1.5±0.2 


283 


R121 


9 h 32 m 10f03 


+21°30W40 


1.9±0.1 


219 


R7f 


9 h 32 m 10?35 


+21°30'09'.'03 


4.0±0.5 


270 


R12m 


9 h 32 m 10f06 


+21°30W24 


3.8±0.3 


212 


R7g 


9 h 32 m 10?29 


+21°30'09'.'64 


2.5±0.2 


259 


R12n 


9 h 32 m 10 s 01 


+21°30W35 


1.8±0.1 


209 


R7h 


9 h 32 m 10?29 


+21°30'09'.'76 


2.2±0.2 


252 


R12o 


9 h 32 m 10?01 


+21°30W36 


2.6±0.2 


207 


R7i 


9 h 32 m 10?37 


+21°30'10'.'54 


2.3±0.1 


248 


R12p 


9 h 32 m 10?01 


+21°30W23 


2.8±0.2 


204 


R7j 


9 h 32 m 10!31 


+21°30'10'.'24 


2.6±0.2 


246 


R12q 


9 h 32 m 10?12 


+21°30W23 


2.3±0.2 


203 


R7k 


9 h 32 m 10?30 


+21°30'09'.'80 


2.6±0.3 


234 


R12r 


9 h 32 m 10 s ll 


+21°30W89 


2.7±0.4 


198 


R71 


9 h 32 m 10!34 


+21°30'09'.'21 


3.2±0.5 


231 


R12s 


9 h 32 m 10?05 


+21°30W17 


2.1±0.1 


193 


R7m 


9 h 32 m 10?26 


+21°30'10'.'19 


2.8±0.2 


230 


R12t 


9 h 32 m 10?16 


+21°30'05^'67 


3.5±0.4 


187 


R7n 


9 h 32 m 10!35 


+21°30'10'.'06 


2.1±0.3 


227 


R12u 


9 h 32 m 10?15 


+21°30W81 


2.1±0.2 


180 


R7o 


9 h 32 m 10?33 


+21°30'09'.'47 


2.1±0.2 


226 


R12v 


9 h 32 m 10 s 14 


+21°30W07 


1.7±0.2 


179 


R7p 


9 h 32 m 10?33 


+21°30'10'.'25 


1.9±0.1 


225 


R4 


9 h 32 m 10?21 


+21°29'59y66 


3.4±0.1 


3772 


R7q 


9 h 32 m 10?32 


+21°30'10'.'56 


1.7±0.1 


222 


R4a 


9 h 32 m 10f08 


+21°30W41 


1.9±0.2 


727 


R7r 


9 h 32 m 10!32 


+21°30'09'.'35 


1.9±0.2 


217 


R4b 


9 h 32 m 10?15 


+21° 29' 59^23 


1.7±0.2 


711 


R7s 


9 h 32 m 10?36 


+21°30'09'.'27 


3.8±0.5 


215 


R4c 


9 h 32 m 10 s 12 


+21°29'59^'45 


1.9±0.1 


703 


R7t 


9 h 32 m 10f34 


+21°30'10'.'35 


2.1±0.4 


212 


R4d 


9 h 32 m 10?ll 


+21°29'59y43 


2.7±0.2 


673 


R7u 


9 h 32 m 10?38 


+21°30'10'.'57 


2.9±0.4 


209 


R4e 


9 h 32 m 10?21 


+21°30W93 


3.4±0.4 


476 


R7v 


9 h 32 m 10?39 


+21°30'10'.'60 


2.2±0.2 


208 


R4f 


9 h 32 m 10?14 


+21°29'59y88 


2.9±0.1 


429 


R7w 


9 h 32 m 10?34 


+21°30'08'.'86 


1.9±0.2 


208 


R4g 


9 h 32 m 10 s 17 


+21°30W05 


2.9±0.4 


310 


R7x 


9 h 32 m 10?32 


+21°30'10'.'04 


2.4±0.2 


207 


R4h 


9 h 32 m 10?16 


+21°30W35 


1.5±0.2 


279 


R7y 


9 h 32 m lK34 


+21°30'08'.'95 


1.7±0.1 


206 


R4i 


9 h 32 m 10?21 


+21°30W09 


2.0±0.3 


274 


R7z 


9 h 32 m 10!28 


+21°30'10'.'30 


2.2±0.1 


204 


R4j 


9 h 32 m 10?18 


+21° 29' 59^54 


2.4±0.2 


271 


N 


9 h 32 m 10?14 


+21°30'03'.'04 


3.8±0.1 


870 


R4k 


9 h 32 m 10 s 27 


+21°29'58y99 


2.0±0.2 


253 


Na 


9 h 32 m 10!ll 


+21°30'02'.'88 


1.5±0.2 


186 


R41 


9 h 32 m 10?10 


+21° 29' 59^78 


2.1±0.3 


246 


Nb 


9 h 32 m 10 s 12 


+21°30'02'.'75 


2.1±0.2 


181 


R4m 


9 h 32 m 10?14 


+21°29'59y52 


2.2±0.3 


245 













the different assumed distance. We also corrected for inter- 
nal extinction using the colour excess [E(B-V)] estimated 
by |Perez^Q lca (1996) from optical spectroscopy and assum- 
ing the galactic extinction law of Miller & Mathews ( 1972 1 



with 7?„=3.2. Planesas et al. (1997) estimated a diameter 
of 2.0arcsec for regions Rl, R2, R7 and the nucleus and 
2.4 arcsec for region R4. In the case of R1+R2 we added their 
Ha luminosities. No values are found in the literature for the 
Ha luminosity of region X. Our derived values of Q(Ho) con- 
stitute lower limits since we have not taken into account the 
fraction of photons that may have been absorbed by dust or 
may have escaped the region. 

Once calculated the number of Lyman continuum pho- 
tons, the masses of the ionizing star clusters, Mi on , have 
been derived using the solar metallicity single burst models 



ionizing photons per unit mass, [Q{Hq) / Mi 0n \. A Salpeter 
initial mass function ( Salpeter|1955 IMF) has been assumed 
with lower and upper mass limits of 0.8 and 120 Mq. In or- 
der to take into account the evolution of the Hn region, we 
have made use of the fact that a relation exists between 
the degree of evolution of the cluster, as represented by the 
equivalent width of the H/3 emission line, and the number of 
Lyman continuum photons per unit solar mass (e.g. |DIaz 



|et al.|2000[ ). We have measured the EW(H/3) from our spec- 



tra (see Tab le [6]) f ollowing the same procedure as in Hagele 



et al. (2006 20081, that is defining a pseudo-continuum to 



by Garcia- Vargas et al. (1995) which provide the number of 



take into account the absorption from the underlying stellar 
population. This procedure in fact may underestimate the 
value of the equivalent width, since it includes the contribu- 
tion to the continuum by the older stellar population (see 
discussions in Diaz et al. 2007 and Dors et al. 2008). The 
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Table 5. Dynamical masses. 



Region 


Slit 


M» 


error(%) 


Region 


Slit 


M* 


error (%) 


Region 


Slit 


M* 


error (%) 


Rla 


SI 


40±5 


12 


R12b 


S2 


71±6 


9 


R4s 


S2 


36±6 


17 


Rib 


SI 


39±6 


17 


R12c 


S2 


71 ±6 


9 


R4t 


S2 


24±4 


17 


Rlc 


SI 


37±5 


13 


R12d 


S2 


59±6 


9 


R4u 


S2 


33±6 


18 


Rid 


SI 


52±6 


12 


R12e 


S2 


45±5 


10 


R4v 


S2 


29±5 


19 


R2 


SI 


46±6 


12 


R12f 


S2 


51±5 


10 


R4w 


S2 


27±4 


16 


R2a 


SI 


47±6 


12 


R12g 


S2 


89±9 


10 


R4x 


S2 


43±8 


18 


R2b 


SI 


47±7 


15 


R12h 


S2 


101±13 


13 


R4y 


S2 


26±5 


20 


R2c 


SI 


46±6 


12 


R12i 


S2 


89±11 


12 


R4z 


S2 


27±5 


20 


R2d 


SI 


44±7 


16 


R12j 


S2 


91±13 


14 


R4a 


S2 


27±5 


20 


R12a 


SI 


47±6 


12 


R12k 


S2 


113±16 


14 


R4sum 


S2 


853±28 


3 


R12b 


SI 


63±7 


11 


R121 


S2 


53±5 


10 


R7 


SI 


26±5 


18 


R12c 


SI 


63±7 


11 


R12m 


S2 


107±12 


11 


R7a 


SI 


26±5 


19 


R12d 


SI 


52±6 


12 


R12n 


S2 


50±5 


10 


R7b 


SI 


32±6 


17 


R12e 


SI 


40±5 


12 


R12o 


S2 


72±8 


11 


R7c 


SI 


26±5 


20 


R12f 


SI 


45±5 


12 


R12p 


S2 


79±8 


10 


R7d 


SI 


25±5 


19 


R12g 


SI 


79±10 


12 


R12q 


S2 


65±7 


11 


R7e 


SI 


14±3 


22 


R12h 


SI 


89±14 


15 


R12r 


S2 


77±12 


15 


R7f 


SI 


38±8 


21 


R12i 


SI 


79±11 


14 


R12s 


S2 


59±6 


9 


R7g 


SI 


24±5 


19 


R12j 


SI 


81±13 


16 


R12t 


S2 


100±13 


13 


R7h 


SI 


21±4 


19 


R12k 


SI 


100±16 


16 


R12u 


S2 


59±7 


12 


R7i 


SI 


22±4 


18 


R121 


SI 


47±6 


12 


R12v 


S2 


47±6 


13 


R7j 


SI 


25±5 


19 


R12m 


SI 


94±13 


13 


R12sum 


S2 


2054±45 


2 


R7k 


SI 


25±5 


21 


R12n 


SI 


44±5 


12 


R12 (adopted) 




1935±66 


3 


R71 


SI 


31±7 


23 


R12o 


SI 


64±8 


13 


R4 


S2 


46±6 


13 


R7m 


SI 


27±5 


19 


R12p 


SI 


70±9 


13 


R4a 


S2 


26±4 


16 


R7n 


SI 


20±5 


22 


R12q 


SI 


58±8 


14 


R4b 


S2 


23±4 


17 


R7o 


SI 


20±4 


20 


R12r 


SI 


68±12 


18 


R4c 


S2 


26±4 


14 


R7p 


SI 


19±3 


18 


R12s 


SI 


52±6 


12 


R4d 


S2 


37±5 


15 


R7q 


SI 


17±3 


18 


R12t 


SI 


88±14 


16 


R4e 


S2 


46±8 


17 


R7r 


SI 


18±4 


20 


R12u 


SI 


52±7 


14 


R4f 


S2 


40±5 


13 


R7s 


SI 


36±8 


22 


R12v 


SI 


42±7 


16 


R4g 


S2 


40±7 


19 


R7t 


SI 


21±5 


25 


R12sum 


SI 


1816±48 


3 


R4h 


S2 


20±4 


19 


R7u 


SI 


28±6 


22 


Rla 


S2 


45±5 


10 


R4i 


S2 


27±5 


20 


R7v 


SI 


21±4 


19 


Rib 


S2 


44±6 


14 


R4j 


S2 


32±5 


15 


R7w 


SI 


19±4 


20 


Rlc 


S2 


41±4 


10 


R4k 


S2 


27±4 


16 


R7x 


SI 


23±4 


19 


Rid 


S2 


59±6 


9 


R41 


S2 


29±5 


19 


R7y 


SI 


17±3 


18 


R2 


S2 


52±5 


10 


R4m 


S2 


30±6 


18 


R7z 


SI 


21±4 


18 


R2a 


S2 


53±5 


10 


R4n 


S2 


20±5 


24 


R7sum 


SI 


642±26 


4 


R2b 


S2 


53±7 


13 


R4o 


S2 


31±6 


18 


N 


S2 


112±10 


9 


R2c 


S2 


52±5 


10 


R4p 


S2 


22±4 


18 


Na 


S2 


43±7 


16 


R2d 


S2 


50±6 


13 


R4q 


S2 


36±6 


17 


Nb 


S2 


62±8 


13 


R12a 


S2 


53±5 


10 


R4r 


S2 


23±4 


17 


Nsum 


S2 


217±15 


7 



masses in 10 b M@. 



derived masses for the ionizing population of the observed 
CNSFRs are given in column 8 of Table [6] and are between 
1 and 4 per cent of the dynamical mass (see column 11 of 
the table). 



The amount of ionized gas (Mhii) associated to each 
star-forming region complex can also be derived from the 
Ha luminosities using the electron density (N e ) dependency 
relation given by |Macchetto et aT. ( 1990 1 for an electron 
temperature of 10 K. The electron density for each region 
(obtained from the [Sn] AA 6717 / 6731 A line ratio) has been 
taken from| Di'az et al.| ( |2007| ) for R1+R2 and R4 and |Perez- 



Qlea| (1996 
300 cm^ 
for the nucleus 



for R7 (see Table [6| . A value of N e equal to 
typical of nuclear Hn regions, has been assumed 



6 DISCUSSION 

6.1 Star and gas kinematics 

The relations between the velocity dispersions of g 
measured from the H/3 (upper panel) and the [Oiii] (lower 
panel) emission lines, and stars, as measured from the IR 
Can triplet, are shown in Fig. [9] In both panels the straight 
line shows the one-to-one relation. In the upper panel, red[^] 
circles show the gas velocity dispersion measured from the 
H/3 emission line using a single Gaussian fit. Orange squares 
and magenta downward triangles show the values measured 
from the broad and narrow components respectively using 
two-component Gaussian fits. The deviant point, marked 
with arrows, corresponds to the galaxy nucleus, whose spec- 
tra have a low signal-to-noise ratio and for which the fits 



3 In all figures, colours can be seen in the electronic version of 
the paper 
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Table 6. Physical parameters. 



Region 


IWHa) a 


E(B-V) 6 


c(Ha) 


L(Ha) 


Q(H ) 


EW(H/3) 


M ion 


N e 


Mhii 


M icm /M« 
(per cent) 


R1+R2 


20.3 


0.53 


0.52 


66.3 


48.7 


12.1 


18.9 


280 c 


0.79 


1.0 


R4 


8.9 


0.66 


0.64 


38.9 


28.6 


4.8 


24.6 


270 c 


0.48 


2.9 


K7 


6.4 


0.71 


0.69 


31.3 


23.0 


3.9 


23.6 


350^ 


0.30 


3.7 


N d 


2.8 






2.8 


2.0 


3.8 


2.1 




0.03 


1.0 



Note. Luminosities in 10 38 erg s 1 , masses in 10 5 Mq, ionizing photons in 10 50 photons 1 and densities in cm 
"From Planesas et al. (1997) corrected for the different adopted distance. 
'From Perez-Olea (1996 ). 



Diaz et al. 



1 2007 I 



c From 

d We assume a value of 0.0 and 300 for E(B-V) and N e , respectively. 



do not provide accurate results for the broad component. In 
the lower panel, red upward triangles, orange diamonds and 
magenta left triangles correspond to the values obtained by 
a single Gaussian fit, and to the broad and narrow compo- 
nents of the two Gaussian fits, respectively. 

In general, the H/3 velocity dispersions of the CNS- 
FRs of NGC 2903 derived by a single Gaussian fit are lower 
than the stellar ones by about 25 kms -1 , except for regions 
R7 and X, for which these values are very similar. These 
two regions also have the lowest velocity dispersions and in 
the first of them the values derived using the two differ- 
ent fitting procedures are very similar. On the other hand, 
a much better agreement between the velocity dispersions 
of stars and the broad component of H/3 for the CNSFRs 
is found. The narrow component shows velocity dispersions 
even lower than those obtained by single Gaussian fits, and 
similar to each other in all cases, with an average value equal 
to 24.5 ± 2.2 km s , with the error given by the dispersion 
of the individual values. This average is represented as a 
magenta horizontal line in the upper panel of Fig. [9] with 
its error as magenta dashed lines. 

The [Oiii] lines however, show a behaviour differ- 
ent from the H/3 lines. In this case, for all the CNS- 
FRs, the gas velocity dispersions derived by a single Gaus- 
sian fit are higher by different amounts, between 7 and 
32 kms -1 . When a two-component fitting procedure is 
used, the narrow component shows the same average value 
found for the H/3 line, although with a larger dispersion 
(a g — 26.2 ± 4.2 kms" 1 ), while the broad component shows 
values even larger than those derived from single Gaussian 
fits. It is therefore apparent that, if there are two kinemati- 
cally different components in the gas, the narrow one seems 
to dominate the H/3 line while the broad one seems to dom- 
inate the [Oiii] ones. In fact, the ratio between the fluxes in 
the narrow and broad components is between 0.65 and 0.95 
for the H/3 line (except for the weakest knot R7, for which 
it is about 1.46) and decreases to between 0.06 and 0.24 for 
the [Oiii] line. 

These results are similar to those found for the previ- 
ously analysed spiral galaxy NGC 3351 (see Paper I). The 
single Gaussian fit yields velocity dispersions which are lower 
than measured for the stars in the case of H/3 and higher in 
the case of [Oiii] by about 20 kms - . 

When two widths Gaussian fits are used, the narrow 
components show a relatively constant value. The average 
value of this velocity for the CNSFRs analysed in both galax- 
ies, NGC 3351 and NGC 2903, amounts to 23.2± 1.7kms -1 
and 25.1±2.0kms -1 for the H/3 and the [Oiii] lines respec- 



tively. Also, for the CNSFRs in the two galaxies, velocity 
dispersions associated with the broad component of H/3 are 
in very good agreement with the stellar ones, except for one 
deviant region in NGC 3351 (R5). For the [Oiii] emission, 
the broad component shows velocities larger than the stel- 
lar ones by up to about 30 kms -1 . 

If the narrow component is identified with ionized gas 
in a rotating disc, therefore supported by rotation, then the 
broad component could, in principle, correspond to the gas 
response to the gravitational potential of the stellar cluster, 
supported by dynamical pressure, explaining the coincidence 
with the stellar velocity dispersion in the case of the H/3 line 
(see |Pizzella et al.|2004| and references therein). The veloc- 
ity excess shown by the broad component of [Oiii] could be 
identified with peculiar velocities in the high ionization gas 
related to massive star winds or even supernova remnants. 

While two velocity components are clearly identified in 
the H/3 line, the weakness of the [Oiii] lines renders this 
identification more uncertain. To test the possibility of find- 
ing a spurious result due to the low signal-to-noise ratio in 
the weak [Oiii] emission lines, we generated a synthetic spec- 
trum with the measured characteristics of radial velocity and 
FWHM adding an artificial noise with a rms twice the ob- 
served one. We then measured the synthetic spectra using 
the same technique used for the real data, that is a dou- 
ble Gaussian fitting with the same initial parameters. In all 
cases we obtained the same result within the observational 
errors. 

Our interpretation of the emission line structures, in the 
present study and in Paper I, parallels that of the studies 
of Westmoquette and collaborators (see for example West- 
moquette et al. 2007a,b). They observed a narrow (~35- 
100 kms -1 ) and a broad (~ 100-400 kms -1 ) component to 
the Ha line across all their four fields in the dwarf galaxy 
NGC 1569. They conclude that the most likely explanation 
of the narrow component is that it represents the general dis- 
turbed optically emitting ionized interstellar medium (ISM), 
arising through a convolution of the stirring effects of the 
starburst and gravitational virial motions. They also con- 
clude that the broad component results from the highly tur- 
bulent velocity field associated with the interaction of the 
hot phase of the ISM (material that is photo-evaporated or 
thermally evaporated through the action of the strong am- 
bient radiation field, or mechanically ablated by the impact 
of fast-flowing cluster winds) with cooler gas knots, setting 



up turbulent mixing layers (e.g. Begelman & Fabian 1990 
|Slavin et al.|1993| . However, our broad component velocity 
dispersion values derived from the H/3 emission lines resem- 
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Figure 9. Upper panel: relation between velocity dispersions of 
the gas (derived from H/3) and stars (CaT) for the CNSFRs and 
the nucleus of NGC 2903. Symbols are as follows: single Gaus- 
sian fit, (red) circles; two Gaussian fit, broad component, (or- 
ange) squares; narrow component, (magenta) downward trian- 
gles. Lower panel: as the upper panel for the [Om] line. (Red) 
upward triangles correspond to the estimates using a single Gaus- 
sian fit, (orange) diamonds represent the broad components of the 
two Gaussian fit and (magenta) left triangles, the narrow com- 
ponents. The (magenta) line represents the average velocity dis- 
persion of the narrow component of the gas (H/3 upper and [Om] 
lower panel) for the CNSFRs, and the (magenta) dashed lines 
represent their estimated errors. [See the electronic edition of the 
Journal for a colour version of this figure.] 



ble more their narrow component values, since in all the 
CNSFRs the former are significantly lower than 100 km s -1 
and show similar values to the stellar velocity dispersions. 
The [Om] emission line on the other hand, shows a be- 
haviour akin to that described by Westmoquette and col- 
laborators. 

There are other studies that identified an underlying 



broad component to the recombination emission lines such 
as iDfaz et all (119871); iTerlevich et all (119961) in the M33 



giant Hll region NGC 604; Chu fc Kennicutt Jr. (19941 



Melnick, Tenorio-Tagle & Terlevich (19991 in the central 



region of 30 Doradus; Mendez & Esteban (19971 in four 



Wolf-Rayet galaxies, and Homcicr & Gallagher ( 1999 ) in the 



starburst galaxy NGC 7673. More recently, Westmoquette et 
al. (2007c, 2009) found a broad feat ure in the Ha emiss ion 



Ostlin et al. 



(2007) in 



lines in the starburst core of M82 
the blue compact g alaxy ESQ 338-IG04; |Sidoli et al.| ( |2006 l 
Firpo et all (120091) in giant extragalactic Hll regions, and 



Hagele| ( |2008| ); |Hagele et al.| ( |2009| ) in CNSFRs of early type 



spiral galaxies. The first two and the last four studies also 
found this broad component in the forbidden emission lines. 

Single Gaussian fits yield the same velocity dispersion 
in the nucleus for both emission lines, H/3 and [Om] 5007 A, 
and close to the stellar value. If a two-component fit is 
used, the low singnal-to-noise ratio in the H/3 line, com- 
bined with the underlying absorption, yields very uncertain 
values for the broad component and a value for the nar- 
row component which is closer to the stellar one than in the 
case of the CNSFRs. The results for the [Om]5007A line 
are similar to those found for the CNSFRs. Fig. [To] shows 
the single and two-component Gaussian fitting for the H/3 
and [Om] 5007 A lines. The H[3 line shows a blue assymetry 
which might correspond to a low intensity broad component. 
This component is barely detectable (if at all present) in the 
[Om] 5007 A line. 

If this decomposition of the emission lines is correct, 
the profiles of the nuclear lines show a behaviour similar 
to that described by Westmoquette and collaborators, with 
the velocity dispersion value of the broad component of H/3 
> 100 km s" 1 



Tschoke et al. 



(2003) 



Using ROSAT data, 
found a very soft X-ray emission feature that extends to 
2' (~ 5 kpc) from the nucleus to the west. They concluded 
that the existence of a galactic wind and its interaction with 
the surrounding intergalactic medium is the most plausi- 
ble source for this soft emission. The H/6eia broad feature 
thus, could be consistent with the presence of a nuclear re- 
gion galactic wind, which could arise from gas that is mixing 
into a high sound speed, hot gas phase, whose presence dom- 
inates X-ray images, as in NGC 1569 (Westmoquette et al. 
2007a,b). 

The presence of two different gaseous components in 
our analysed CNSFR could have an important effect on the 
classification of the activity in the central regions of galax- 
ies through diagnostic diagrams and the magnitude of this 
effect would increase with decreasing spatial resolution. In 
fact, the [Om]/H/3 ratio is much smaller for the narrow com- 
ponent than for the broad one, and single Gaussian fittings 
provide intermediate values (see Table This seems to 
point to a lower excitation for the lower velocity kinemat- 
ical component. The effect is rather dramatic even for the 
nucleus of NGC 2903 for which the logarithmic [Om]/H/3 
ratio for the low and high velocity components vary from - 
0.91 to -0.22. Unfortunately, we do not have high resolution 
data in the [Nil] - Ha wavelength range and hence we can 
not investigate the complete effect in the excitation plot (see 
Paper I). 

The radial velocities along the slit for each angular po- 
sition of NGC 2903 as derived from the ionized gas emission 
lines, H/3 and [Om] 5007 A, and the stellar CaT absorptions 
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Figure 10. Single and two-component Gaussian fits for the H/3 (left panel) and [Oin] A 5007 A (right panel) emission lines in the nucleus 
of NGC 2903. Thick dashed (red) line and two-component fits as resulting from the IRAF task ngaussfit are superposed. For the two- 
component fit, the dashed-dotted line corresponds to the broad component, the dotted line to the narrow component and the dashed 
line to the sum of both. 



are shown in the upper and lower panels of Fig. [TT] respec- 
tively. The rotation curves seem to have the turnover points 
at the same positions as the star-forming ring, specially for 
the S2 slit position across the nucleus, as found in other 
galaxies (see Telesco & Dccher 1988; Diaz et al.||1999 and 
references therein). For the systemic velocity of NGC 2903, 
the derived values are consistent with those previously ob- 



tained by Planesas et al. (1997) and Haynes et al. (19981, 



and with the velocity distribution expected for this type of 
galaxies (Binney & Tremainc 19871. 

The radial velocities derived from both the stellar CaT 
and the [Oiii]5007A emission line are in good agreement. 
The H/3 line velocities, however, differ from the stellar ones, 
more appreciably in the Norh-West region, in amounts sim- 
ilar to the differences shown by the radial velocities corre- 
sponding to the two Gaussian components (Av n ;,; see Table 
[5]). This is compatible with the assumption of the existence 
of two kinematically different components, the broad one 
dominating the single fit in [Oin] and the narrow one being 
the dominant one in the case of H/3. 

Fig. [12] shows the run of velocity dispersions along 
the slit versus pixel number for slit positions SI and S2 
of NGC 2903, respectively. These velocity dispersions have 
been derived from the gas emission lines, H/3 and [Oin], and 
the stellar absorption ones using 2 and 3px apertures for 
SI and S2, respectively. We have marked the location of 
the studied CNSFR and the galaxy nucleus. We have also 
plotted the stellar velocity dispersion derived for each region 
and the nucleus using the 5 px apertures. The values derived 
with wider apertures are approximately the average of the 
velocity dispersions estimated using the narrower apertures. 
Then, the increase of the apparent velocity dispersions of 
the regions lying on relatively steep parts of the rotation 
curve by spatial "beam smearing" of the rotational velocity 
gradient due to the finite angular resolution of the spectra 
is no so critical for the studied CNSFRs. This is due to the 
turnover points of the rotation curves and the star-forming 
ring seem to be located at the same positions, as was pointed 
out above. The behaviour of the velocity dispersion along 
slit position SI (upper panel) is very complex. The H/3 line 
shows a relatively constant value of the velocity dispersion. 
To the SW of the R1+R2 star forming complex, there is a 



clear systematic difference between the velocities measured 
from the HI recombination line and the stellar absorption 
lines. This systematic difference is significative lower in the 
region where active star formation is taking place. Along 
slit position S2, crossing the galaxy nucleus, we detect rel- 
ative minima in the stellar and gas velocity dispersions at 
the positions of CNSFR R4 and R1+R2. The stellar velocity 
dispersions show also a relative minimum at the location of 
the nucleus although in this case it is accompanied by a rise 
in the velocity dispersion measured from the H/3 line. 



6.2 Star cluster masses 

Unlike the case of the CNSFRs in NGC 3351 studied in Pa- 
per I for which we found that two of the observed regions, R4 
and R5, seem to possess just one knot showing up in the con- 
tinuum image and coincident with the Ha emission, the four 
CNSFRs observed in NGC 2903: R1+R2, R4 and R7 show 
a complex structure at the HST resolution (Fig. [8]l , with a 
good number of subclusters with linear diameters between 
3 and 8 pc. For these individual clusters, the derived upper 
limits to the masses are in the range between 1.4 x 10 6 and 
1.13 x 10 7 M Q (see table |||, with fractional errors between 
about 9 and 25 per cent. The upper limits to the dynam- 
ical masses estimated for the whole CNSFRs ("sum") are 
between 6.4 x 10 7 and 1.9 x 10 8 Mq, with fractional errors 
between about 3 and 4 per cent. The upper limit to the dy- 
namical mass derived for the nuclear region inside the inner 
3.8pc is 1.1 x 10 7 Mq, with a fractional error of about 9 per 
cent. 

The masses of the ionizing stellar clusters of the cir- 
cumnuclear complexes, have been derived from their Ha lu- 
minosities under the assumption that the regions are ion- 
ization bound and without taking into account any pho- 
ton absorption by dust. A range of masses between 1.9 and 
2.5 x 10 6 M© for the star- forming regions, and 2.1 x 10 5 M Q 
for the nucleus are found (see Table [6|. In column 11 of 
Table [6] we show a comparison (in percentage) between the 
ionizing stellar masses of the circumnuclear regions and their 
dynamical masses. These values are approximately between 
1-4 per cent for the CNSFRs, and 1 per cent for the nucleus 
of NGC 2903. 
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Figure 11. Radial velocities along the slit versus pixel number 
for each slit position of NGC 2903 (upper panel: SI; lower panel: 
S2) as derived from the gas emission lines (red circles: H/3; upward 
triangles: [Olll]) and the stellar absorption ones (blue squares). 
The individual CNSFRs and the nucleus, "N" , or the closest po- 
sition to it, "d" , are marked in the plots. The dashed-dotted line 
is the systemic velocity of NGC 2903 derived by Planesas et al. 
(1997). The distance in arcsec from the nucleus is displayed in 
the upper x-axis of each panel. [See the electronic edition of the 
Journal for a colour version of this figure.] 
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Figure 12. Velocity dispersions along the slit versus pixel num- 
ber for each slit position of NGC 2903 (upper panel: SI; lower 
panel: S2) as derived from the gas emission lines (small red cir- 
cles: H/3; small yellow triangles: [Olll]) and the stellar absorption 
ones (small blue squares). The stellar velocity dispersions derived 
for each region and the nucleus using the 5 px aperture are also 
plotted with violet squares. The individual CNSFRs and the nu- 
cleus, "N", or the closest position to it, "d", are marked in the 
plots. The distance in arcsec from the nucleus is displayed in the 
upper x-axis of each panel. [See the electronic edition of the Jour- 
nal for a colour version of this figure.] 



Finally, the masses of the ionized gas, also derived from 
their Ha luminosities, range between 3 x 10 4 and 8 x 10 5 M Q 
for the CNSFRs, and 3 x 10 3 M for the nucleus (see Table 
[6). They make up a very small fraction of the total mass 
of the regions. Both the masses of the ionizing stellar clus- 
ters and those of the ionized gas are comparable to those 
derived by Gonzalez-Delgado et al. ( 1995 1 for the circumnu- 



clear region A in NGC 7714. It should be taken into account 
that the latter have been derived from the Hcv luminosity of 
the CNSFRs assuming a single kinematical component for 
the emission line. If we consider only the broad component 
whose kinematics follows that of the stars in the regions, all 
derived quantities would be smaller by a factor of 2. 

The circumnuclear star formation in NGC 2903 has 



been studied by Alonso-Herrero et al. (20011 using high 



resolution near-IR photometry in the H band and ground 
based IR spectroscopy. Although, in general, the spatial 
distribution of the stellar clusters detected in the IR does 



not coincide with the maxima of the Ho emission, some 
CNSFR show up prominently in both spectral ranges. This 
is the case for our regions R4, R7 and R1+R2 which can 
be identified with regions H8, H4 and H1+H2, respectively, 
in Alonso-Herrero et al. (20011. These authors have esti- 



mated the masses of the young clusters present in the cen- 
tral region of NGC 2903 from the NIR colours using stellar 
population synthesis models corresponding to single bursts 
of star formation with Gaussian FWHM values of 1, 5 and 
100 Myr and a Salpeter IMF with lower and upper mass 
cut-offs of 1 and 80 M Q . Values between 2.1 x 10 8 M and 
3.6 x 1O 8 M are obtained for the central ~625pc region 
without including the young stellar population responsible 
for the bright Hn region emission. This is consistent with 
our upper limit to the mass corresponding to the sum of 
the dynamical mass found for the CNSFR and the nuclear 
clusters (3.6 x 10 s M ). 
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The values of dynamical masses of the individual clus- 
ters derived by us for the CNSFRs of both NGC 2903 and 
NGC 3351 (see Paper I) are larger than those estimated for 
the individual clusters in the nuclear ring of NGC 4314 from 



in the simulations by Bosch & Meza (2001 1 are set for their 



broad band U,B,V,I and Ha photometry with HST (Bene- 



dict et al. 2002 \ by more that an order of magnitude. No 



comparable data exist for the CNSFRs studied here to al- 
low the estimate of photometric M/L ratios. However, the 
fact that the H/3 luminosities of the NGC 4314 clusters are 
lower than those of NGC 2903 by a factor of about 40 gives 
some support to our larger computed dynamical masses. 
These masses are, on average, of the order of the largest 



kinematically derived mass for SSC in M82 (McCrady & 



Graham 



2007 1 from near IR spectroscopy (40 x 10 b M ) but 



also larger (by factors between 4.2 and 34) than the mass 
derived for the SSC A in NGC 1569 by |Ho fc Filippenko] 
1 1996a[ ) from stellar velocity dispersion measurements using 
red (~6000A) spectra (3.3±0.5 x 1O 5 M ). 

It should be recalled that we have estimated the dynam- 
ical masses through the virial theorem under the assumption 
that the systems are spherically symmetric, gravitationally 
bound and have isotropic velocity distribution. We have used 
the stellar velocity dispersions derived from the CaT absorp- 
tions features and the cluster sizes measured from the high 
spatial resolution WFPC2-PC1 HST image. 

Therefore, while the average radius of an individual 
cluster is of the order of 2.8 pc, the stellar velocity dispersion 
corresponds to a much larger region, typically ~ 42 x 76 pc 2 
(1.0 x 1.8 arcsec 2 ) containing several knots. The use of these 
wider size scale velocity dispersion measurements to esti- 
mate the mass of each knot, can lead to an overestimate 
of the mass of the individual clusters, and hence of each 
CNSFR (see Paper I). 

However, as can be seen in the HST-NICMOS image 
(right-hand panel of Fig. []}, the CNSFRs of NGC 2903 
show up very prominently in the near-IR and dominate the 
light inside the apertures observed (see|Alonso-Herrero et al!] 



2001 for a detailed analysis of the IR images). Therefore the 



assumption that the light at the CaT wavelength region is 
dominated by the stars in the clusters seems to be justi- 
fied. The IR CaT is very strong, in fact the strongest stellar 
feature in very young clusters, i.e. older than 4Myr (Ter- 



|levich et al.|[T990| ). Besides, we detect a minimum in the 
velocity dispersion at the position of the clusters (see Fig. 
121. Nevertheless, we can not assert that we are actually 
measuring their velocity dispersion and thus prefer to say 
that our measurements of cr*, and hence dynamical masses, 
constitute upper limits. Although we are well aware of the 
difficulties, yet we are confident that these upper limits are 
valid and important for comparison with the gas kinematic 
measurements. 

Another important fact that can affect the estimated 
dynamical masses is the presence of binaries among the red 
supergiant and red giant populations from which we have de- 
rived the stellar velocity dispersions. In a recent work, Bosch 
et al. (2009) using GEMINI-GMOS data, have investigated 
the presence of binary stars within the ionizing cluster of 
30 Doradus. From a seven epochs observing campaign they 
have detected a rate of binary system candidates within their 
OB star sample of ~ 50 %. Interestingly enough, this detec- 
tion rate is consistent with a spectroscopic population of 
100 % binaries, when the observational parameters described 



observations. From their final sample of 'single' stars they 
estimated a radial velocity dispersion of 8.3kms _1 . When 
they derived <r» from a single epoch, they found values as 
high as SOkms" 1 , consistent with the values derived from 



single epoch NTT observations by Bosch et al. ( 2001 1 



Although the environment of our CNSFRs is very dif- 
ferent from that of 30 Dor and the stellar components of 
the binary systems studied by Bosch et al. (2009) are very 
different from the stars present in our regions from where 
the CaT arise (red supergiants), this is an illustrative ob- 
servational example of the problem. The orbital motions of 
the stars in binary (multiple) systems produce an overesti- 
mate of the velocity dispersions and hence of the dynamical 
masses. The single-star assumption introduces a systematic 
error that depends on the properties of the star cluster and 
the binary population, with an important effect on the clus- 
ter mass if the typical orbital velocity of a binary component 
is of the order of, or larger than, the velocity dispersion of the 



single/binary stars in the potential of the cluster (Kouwen 



hoven & de Grijs 2008 1. As was pointed out by these authors, 



the relative weights between the single and binary stars in 
the velocity dispersion measurements depend on the binary 
fraction, which, together with the semi-major axis or period 
distribution, are the most important parameters in order 
to determine if the binary population affects the estimated 
dynamical masses. Their simulations indicate that the dy- 
namical mass is overestimated by 70 % , 50 % and 5 % for 
a measured stellar velocity dispersion in the line of sight of 
lkms -1 , 2kms~ 1 and lOkrns -1 respectively. They there- 
fore conclude that most of the known dynamical masses of 
massive star clusters are only mildly affected by the presence 
of binaries. Hence, although the binary fraction of the red 
supergiants and red giants in this type of circumnuclear and 
typically high metal rich environment is unknown, for our 
clusters, where the smallest estimated velocity dispersion is 
37kms _1 , we can assume that the contribution of binaries 
to the stellar velocity dispersions is not important. 



7 SUMMARY AND CONCLUSIONS 

We have measured gas and stellar velocity dispersions in 
four CNSFRs and the nucleus of the barred spiral galaxy 
NGC 2903. The stellar velocity dispersions have been mea- 
sured from the CaT lines at AA8494, 8542, 8662 A, while 
the gas velocity dispersions have been measured by Gaus- 
sian fits to the H/3 A 4861 A and the [Om] A5007A emission 
lines on high dispersion spectra. 

Stellar velocity dispersions are between 37 and 
65kms _1 . These values are about 25kms _1 larger than 
those measured for the ionized gas from the H/3 emission 
line using a single Gaussian fit. On the other hand, single 
Gaussian fits to the [Om] 5007 A line yield velocity disper- 
sions close to the stellar ones, and, in some cases, somewhat 
larger. However, the best Gaussian fits involved two differ- 
ent components for the gas: a "broad component" with a 
velocity dispersion similar to that measured for the stars, 
and a "narrow component" with a velocity dispersion lower 
than the stellar one by about SOkms" 1 . This last velocity 
component shows a relatively constant value for the two gas 
emission lines, close to 24kms _1 for all the studied CNSFRs 
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(including the regions belonging to NGC 3351 from Paper I). 
We find a radial velocity shift between the narrow and the 
broad component of the multi- Gaussian fits in the CNSFRs 
of NGC 2903 that varies between -10 and 35kms~ . 

When plotted in an [Om]/H/3 versus [Nil] /Ha dia- 
gram, the two systems are clearly segregated for the high- 
metallicity regions of NGC 2903, with the narrow compo- 
nent showing lower excitation and being among the lowest 
excitation line ratios detected within the SDSS data set of 
starburst systems. Unfortunately, no comparable informa- 
tion exists for the [Nil] /Ha ratio and hence we can not eval- 
uate any possible contribution by shocks. 

The upper limits to the dynamical masses estimated 
from the stellar velocity dispersion using the virial theo- 
rem for the CNSFRs of NGC 2903 are in the range between 
6.4 x 10 7 and 1.9 x 10 8 M Q and is 1.1 x 10 7 M for the nu- 
clear region inside 3.8 pc. Masses derived from the H/9 veloc- 
ity dispersions under the assumption of a single component 
for the gas would have been underestimated by factors be- 
tween 2 and 4 approximately. The upper limits to the derived 
masses for the individual clusters are between 1.4 x 10 6 and 
1.1 x 1O 7 M0. These values are between 4.2 and 34 times 
the mass deriv ed for the SSC A in NGC 1569 by |Ho fc Fili"p~J, 
penko ( 1996a I and larger than other kinematically derived 



SSC masses. 

Masses of the ionizing stellar clusters of the CNSFRs 
have been derived from their Ha luminosities under the as- 
sumption that the regions are ionization bound and with- 
out taking into account any photon absorption by dust. 
For the star forming complexes of NGC 2903 these masses 
are between 1.9 and 2.5 x 1O B M , and is 2.1 x 10 5 M for 
the galaxy nucleus. Therefore, the ratio of the ionizing stel- 
lar population to the total dynamical mass is between 0.01 
and 0.04. These values of the masses of the ionizing stellar 
clusters of the CNSFRs are comparable to that derived by 



Gonzalez-Delgado et al. ( 1995 1 for the circumnuclear region 



A in NGC 7714. 

Masses for the ionized gas, also derived from their Ha 
luminosities, vary between 3.0 x 10 4 and 7.9 x 10 5 M for 
the star forming regions and is 3 x 10 3 M for the nucleus of 
NGC 2903. These values are also comparable to that derived 



by Gonzalez-Delgado et al. ( 1995 1 



It is interesting to note that, according to our findings, 
the SSC in CNSFRs seem to contain composite stellar pop- 
ulations. Although the youngest one dominates the UV light 
and is responsible for the gas ionization, it constitutes a few 
per cent of the total mass. This can explain the low equiv- 
alent widths of emission lines measured in these regions. 
This may well apply to other SSC and therefore conclusions 
drawn from fits of single stellar population (SSP) models 
should be taken with caution (e.g. |McCrady et al.||2003| 



Larsen et al. 20041. Furthermore, the composite nature of 



the CNSFRs means that star formation in the rings is a 
process that has taken place over time periods much longer 
than those implied by the properties of the ionized gas. 

The observed stellar and [Oiii] radial velocities of 
NGC 2903 are in good agreement, while the H/3 measure- 
ments show shifts similar to those found between the nar- 
row and the broad components. This different behaviour can 
be explained if the positions of the single Gaussian fits are 
dominated by the broad component in the case of the [Om] 
emission line and by the narrow one in the case of H/9. The 



rotation curve corresponding to the slit position through 
the nucleus shows maximum and minimum values at the 
positions of the circumnuclear regions, as observed in other 
galaxies with CNSFRs. 

The existence of more than one velocity component in 
the ionized gas corresponding to kinematically distinct sys- 
tems, deserves further study. Several results derived from 
the observations of the different emission lines could be af- 
fected, among others: the classification of the activity in the 
central regions of galaxies, the inferences about the nature 
of the source of ionization, the gas abundance determina- 
tions, the number of ionizing photons from a given region 
and any quantity derived from them. To disentangle the ori- 
gin of these two components it will be necessary to map 
these regions with high spectral and spatial resolution and 
much better S/N ratio in particular for the 2+ lines. High 
resolution 3D spectroscopy with IFUs would be the ideal 
tool to approach this issue. 
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